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We discuss the universal nature of relaxation in isolated many-body quantum systems subjected 
to global and strong periodic driving. Our rigorous Floquet analysis shows that the energy of the 
system remains almost constant up to an exponentially long time in frequency for arbitrary initial 
states and that an effective Hamiltonian obtained by a truircatioir of the Floquet-Magnus expansion 
is a quasi-conserved quantity in a long timescale. These two general properties lead to iirtriguiirg 
classification on the initial stage of relaxation, one of which is similar to the prethermalization 
phenomenon in nearly-integrable systems. 


Introduction .— In periodically driven many-body 
quantum systems, excited states as well as the ground 
state participate in the dynamics, and nontrivial macro¬ 
scopic phenomena can appear. Recent years have wit¬ 
nessed remarkable experimental developments, such as 
the discoveries of the Higgs mode in the oscillating or¬ 
der parameter of the superconducting material under a 
Terahertz laser [1], and the Floquet topological states in 
the periodically driven cold atom [2-5]. Periodic driving 
in isolated quantum systems sometimes generates unex¬ 
pected dynamical phenomena even if the instantaneous 
Hamiltonian at each time step is simple. To name only a 
few, dynamical localization [6-8] , coherent destruction of 
tunneling [7-9], dynamical freezing [10, 11], and dynam¬ 
ical phase transitions [12-14] are remarkable far-from- 
equilibrium phenomena that cannot be captured within 
linear-response analysis. 

On the other hand, as recently discussed in the context 
of thermalization, careful consideration is necessary on 
the true steady state in driven many-body systems [15- 
19]. Thermalization in isolated quantum systems has 
become one of critical subjects in modern physics [20- 
25]. The first study was made by von Neumann early 
in 1929 [26], and now we are on the new stage by in¬ 
corporating many concepts including quantum entangle¬ 
ment [27] and experiments [28]. In the case without driv¬ 
ing fields, the notion of the eigenstate thermalization hy¬ 
pothesis (ETH) is a key idea [20, 21, 23, 26] that states 
that each energy eigenstate is indistinguishable from the 
microcanonical ensemble with the same energy. As a 
generalization of ETH to periodically driven systems, the 
Floquet ETH was proposed, which states that all the Flo¬ 
quet eigenstates look the same and are indistinguishable 
from the infinite-temperature (i.e., completely random) 
state [16, 19, 29-31]. This leads to the conclusion that 
in general periodically driven many-body systems will 
eventually reach the steady state of infinite temperature, 
although several exceptions exist [15, 17-19]. 

The question that follows is on the time scale to reach 


the steady state. Recent experiments seem to urge us 
to clarify the general aspects of the time scale especially 
for the strong amplitude of global driving, where non¬ 
trivial transient dynamics is anticipated. We note that 
most nontrivial dynamical phenomena in driven systems 
are far-from-equilibrium effects that cannot be analyzed 
within linear-response analysis. Hence, in this paper, we 
for the hrst time aim to find the universal nature of the 
relaxation to the steady state under strong and global 
driving. This direction is clearly crucial for a deeper un¬ 
derstanding of thermalization and for analyzing the sta¬ 
bility of transient dynamics in experiments. 

For this aim, we focus on the Floquet Hamiltonian Hp 
which plays a central role in periodically driven systems: 

where Hit) is the Hamiltonian of the system, T is the 
time-ordering operator, and T is the period of the driv¬ 
ing {h = 1 throughout this paper). The Floquet Hamil¬ 
tonian is an effective Hamiltonian that contains full in¬ 
formation on the stroboscopic dynamics. The Floquet- 
Magnus (FM) expansion is a formal expression for the 
Floquet Hamiltonian: Hp = 

explicit form of Ulm is given in Eq. (8) below. However, it 
has recently been recognized that using full series expan¬ 
sion is problematic since it is not convergent in general. 
The convergence radius shrinks as the system size in¬ 
creases [33] . Instead we here use the technique of trunca¬ 
tion in the FM expansion, which was recently developed 
for describing the Floquet Hamiltonian for transient time 
scales [34, 35]: 

n 

jjin) ^ ■ ( 2 ) 

m—0 

(n) 

Here, Hp is the nth order truncated Floquet Hamil¬ 
tonian. There are several studies which show that the 
time-evolution by the truncated Floquet Hamiltonian is 
reliable up to a certain long time r for the driving with 
high frequency oj = jT [33]; r ~ for the Friedrichs 
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model on the continuous space [34] and r ^ exp[C>(a;)] 
for lattice systems when driving is local [35] or interac¬ 
tions are short-ranged [35-37]. In this paper, we use the 
truncation technique for high frequency driving. 

With this technique, two findings are mainly presented. 
We show as the first result that in the case of a high- 
frequency driving, the truncated Floquet Hamiltonian is 
a quasi-conserved quantity (a quantity that is almost con¬ 
served in a long timescale). We also show as the second 
result that energy absorption rate per one site is bounded 
for an arbitrary amplitude of driving and for arbitrary 
initial states: 

E/N < {Nv/N) exp[-0(u;/5)], (3) 

where E and N are respectively the total energy and the 
number of lattice sites, and g is the maximum energy 
per one site. The driving field is applied to Ny sites. 
This provides a criterion on stability of transient quan¬ 
tum dynamics in experiments. These two findings lead to 
intriguing classification on the relaxation processes, one 
of which is similar to the prethermalization phenomenon 
seen in non-driven nearly-integrable systems [38-41] , see 
Refs. [42, 43] for recent relevant numerical calculations. 

Setup and numerical example .— We consider a quan¬ 
tum spin system defined on a lattice with N sites in ar¬ 
bitrary dimension, whose Hamiltonian is written as 

H{t) = Ho + Vit ). (4) 

The driving field V{t) is applied to Nv{< N) sites and 
satisfies the periodicity in time V{t) = V{t + T) with 
zero average over the single period. We mainly focus on 
the regime of high-frequency to = 2 tt/T. Each lattice site 
i = 1, 2,... has its own spin. The basic assumption on 
the Hamiltonian is that it is expressed in the form of 

H{t)= (5) 

X:|X|<fe 

where X = {ii,i 2 , ■ ■ ■ ,i\x\} is a set of the lattice sites 
with |Xj being the number of sites in X, and hx{t) is an 
operator acting on the sites in X. In addition, we assume 
that the single site energy is bounded in the sense that 

for any site i, E ll'■x(^)ll<9 (6) 

X-.X^i 

with some fixed positive constant g, where jj • jj denotes 
the operator norm. 

The form of Eq. (5) means that the Hamiltonian con¬ 
tains at most fc-body interactions. For most physical ap¬ 
plications, we can consider the case of fc = 2. In the 
case of spin-(1/2) systems, the most general form of the 
Hamiltonian (5) with k = 2 is 

N N 

77(f) =^B.(t).<T, + ^ ^ J^^t)aTa], (7) 

i—1 idj (y,'i=x,y,z 
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FIG. 1. (color online) Numerical demonstration of 
prethermalization-like phenomenon, (a): Relaxation in the 
long timescale, (b): initial relaxation, (c): transient time- 
evolution after the initial one. Parameters are {J, B^, Bz) = 
(1,0.9045, 0.8090) and N = 24. The dotted line in (b) is the 
expectation value in the equilibrium state of at the in¬ 
verse temperature /3 = 0.85, which is determined from the 
expectation value of at t = 10. 


where cr^ = (erf, cr/, cr/) is the Pauli matrix of ith spin, 
Bi{t) is the local magnetic field at ith site, and [t) de¬ 
notes the interaction between 7th and jth spins. We can 
explicitly confirm that this Hamiltonian can be brought 
into the form of Eq. (5) by putting /ijjj. = Bi ■ (Ji and 

To make clear physical phenomena that we address, 
we show a numerical example with a toy model that has 
been used to show the Floquet ETH in [31]. We consider 
dynamics over one cycle, taking H{t) = Hz for the first 
half period and 77(f) = Hx for the second half period, 
where Hz = ^*=1 + Bzcrf] with the periodic 

boundary condition and Hx = BxJ2f=i’^i- We calculate 
the time-evolution of the z-component of the hrst spin erf 
setting all spin-down state as the initial state. Floquet 
ETH implies that a steady state in the long-time limit is 
a random state, and hence, when it is satisfied, the expec¬ 
tation value of any local spin operator eventually reaches 
zero. In Fig.l, the time-evolution for a sufficiently large 
system size is shown. Figure 1(a), (b) and (c) are re¬ 
spectively time-evolution of erf in the large timescale, 
the initial stage, and the transient timescale. Fig. 1(a) 
shows a vanishing expectation value that is a clear indi¬ 
cation of the Floquet EHT. Crucial observation is that 
after the initial relaxation (Fig.1(b)), the expectation 
value is almost constant for finite transient timescales, 
and the timescales depend on the period T (Fig. 1(c)). 
This implies that the heating process is seemingly sup¬ 
pressed during this timescale. This is somewhat simi- 
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lar to the prethermalization phenomenon in nondriven 
nearly-integrable systems. In experimental situations, 
this transient time behavior is crucial, and hence we ad¬ 
dress the mechanism of the behavior and consider the 
period dependence on the timescale. 


Timescale of the heating process .— We use the FM 
expansion for analyzing the energy absorption and the 
relaxation process. The FM expansion is the formal 
expansion of the Floquet Hamiltonian given hy Hp = 
J2'^=o with Hq = Hq and the nth order coefficient 

for n > 1 being given by [44] 




= E 


i'^{n + 




( 8 ) 


where ct is a permutation and 6[a\ = + 1) ~ 

(j{i)) with O(-) is the step function. It is believed that 
the FM expansion is divergent in many-body interacting 
systems [19, 29, 33]. See the supplementary material for 
the numerical demonstration of the divergence [45] . This 
divergence is not merely a mathematical phenomenon but 
is now thought to be an indication of heating process due 
to periodic driving [19, 29, 33]. 

We define the nth order truncation of the FM expan¬ 
sion as in Eq. (2) and show that for general spin systems 
the timescale of the heating is exponentially slow in fre¬ 
quency. To this end, we start with an intuitive expla¬ 
nation on our analysis. From Eq. (8), has at most 
(n -|- l)fc-spin effective interactions because of the multi¬ 
ple commutators in Eq. (8), which describes the collective 
flip of (n + l)k spins. Since the energy exchange between 
a quantum system and the external periodic field is quan¬ 
tized into integer multiples of uj and the energy of each 
spin is bounded by g, N* ^ uj/g spins must flip coop¬ 
eratively in order to absorb or emit the single “energy 
quantum”. Such a process is taken into account only in 
the terms higher than the noth order in the FM expan¬ 
sion with no ^ N*/k ^ uj/{gk). Indeed, each term of the 
FM expansion is rigorously bounded from above as 


\\^n\\T” < 2gNv 


{2gkTYn\ 
(n -I- 1)2 


( 9 ) 


This is given by estimating norms of the multiple commu¬ 
tators taking account that the Hamiltonian has at most 
fc-body interaction and the energy per site is bounded by 
g [35, 45]. Equation (9) shows that the FM expansion (2) 
looks convergent up to n < no w/ (gk) and 


\\hP - = lVyO(T"+i) (n < no), (10) 


but grows rapidly for n > uq. Therefore, we can elimi¬ 
nate the heating effect most efficiently by truncating the 
FM expansion at n = np. The timescale of the heating is 
thus evaluated by comparing the difference between the 
exact time evolution and the approximate time evolution 
under the npth order truncated Floquet Hamiltonian. It 
is expected that higher-order terms (i.e. simultaneous 


flip of a large number of spins) would matter only in the 
later stage of the time evolution. 

We now make the above argument mathematically rig¬ 
orous. We can prove the following theorem: 

Theorem. The ngth order truncated Floquet Hamilto¬ 
nian is almost conserved up to an exponentially 

long time in frequency in the sense that 

\\H^p°\t) - < 16g^k2-^°Nvt, (11) 

where t = niT with a positive integer m, uq = 
[l/iSgkT) - IJ, and H^p°\t) = W {t)H^p°^Uit) is the 
noth order truncated Floquet Hamiltonian at time t in 
the Heisenberg picture with U{t) = Te~’‘do ). 

This is derived by evaluating the norm of the Dyson- 
expansion for the time-evolution operators in the left 
hand side taking into account that the Hamiltoanin is 
written as Eq. (5) with Eq. (6). See the supplementary 
material for more details on the derivation [45]. Com¬ 
bined with Eq. (10), this theorem leads to 

\\H^p\t)-H^p'>\\ < 16g2fc2-"“7Vyf-hiVyO(r"+i) (12) 

for any n < uq. In particular, by substituting n = 0, we 
obtain 

^\\Ho{t) -Ho\\<^ [I6g^k2--H + 0{T)] . (13) 

It is noted that the term of 0{T) in Eq. (13) is inde¬ 
pendent of t. Thus, the energy density remains con¬ 
stant within a small fluctuation of 0{T) until an ex¬ 
ponentially long time in frequency. Equation (11) pro¬ 
vides a lower bound on the timescale during which H^°^ 
can be an approximately conserved quantity. This quasi- 
conserving property lasts during the timescale larger than 
r ~ 2"“ ~ Similarly, Eq. (13) implies that the 

lower bound of the timescale of heating is an exponen¬ 
tially long time in frequency. This is a main result (3). 
The exponentially long timescale of the energy relaxation 
was shown for short-range interacting spin systems in 
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the linear-response regime in Ref. [46], but it should be 
emphasized that Eq. (13) has been obtained without as¬ 
suming short-range interactions and the linear-response 
argument. See Ref. [47] for a recent numerical result. 

It is remarked that for local driving with Ny < 
a much stronger result was shown in Ref. [35], i.e., 

|| 7 -e-* fo H(s)ds _ II < exp[-0(w)]t (14) 

for t = mT. This inequality implies that for any bounded 
operator that may be highly nonlocal, the FM truncated 
Hamiltonian gives the accurate time evolution up to an 
exponentially long time. In the case of global driving 
Ny oc N, this strong inequality (14) is not satisfied for 
sufficiently large systems, but even in this case, we can 
utilize the finite order truncation of the FM expansion to 
discuss the relaxation process as is argued below. 

Relaxation process. — Our rigorous result enables us 
to discuss possible scenarios on the initial stage of re¬ 
laxation. According to the Floquet ETH [19, 29-31], 
the steady state in the long-time limit induced by the 
Floquet Hamiltonian (1) is a state of infinite tempera¬ 
ture. Full FM series expansion in general diverges in 
large quantum systems and hence it is not useful [33, 35]. 
However, the truncated Floquet Hamiltonian is a 

quasi-conserved quantity and plays a crucial role in the 
relaxation process. 

We make a remark on the degree of nonlocality on 
the quasi-conserved quantity. The noth order truncated 
Floquet Hamiltonian has effective (no -I- l)fc-body inter¬ 
actions, and hence the nonlocality looks large. However, 
for a high frequency driving, higher order contributions 
in the FM expansion are very small, since T is small. The 
dominant contribution is in fact the original Hamiltonian 
Hq. Hence, nonlocality of the truncated Floquet Hamil¬ 
tonian is not very strong. Eigenstates for the truncated 
Floquet Hamiltonian thus should satisfy the usual 

ETH, not the Floquet ETH. In addition, we should note 
that from Eq. (10), ~ for any n < uq. Hence 

these truncated Floquet Hamiltonians are not indepen¬ 
dent but almost the same. Practically one can approxi¬ 
mate the quasi-conserved quantity by Hq{= Hp^). 

Taking account of those, we discuss a scenario on the 
initial stage of relaxation. Since the quasi-conserved 
quantity exists with long lifetime, the system relaxes 
to a quasi-stationary state characterized by the quasi- 
conserved quantity, which will be close to a state corre¬ 
sponding to the (micro) canonical ensemble of the 

effective Hamiltonian H^p°'^ set by the initial state. Ap¬ 
proximately one can use (the equilibrium ensemble 
of Hq) instead of because for any 

n < Hq. 

The initial stage of relaxation can be classified into 
two cases, i.e., (i) the case where the relaxation to the 
quasi-stationary state is faster than the energy relax¬ 


ation, and (ii) the case where both relaxation times are 
comparable. In the case (i), the system first reaches 
the quasi-stationary state, and then relaxes to the true 
steady state. This is highly related to the prethermaliza- 
tion phenomenon in the isolated nearly-integrable sys¬ 
tems [40, 41], where the system first relaxes to a quasi- 
stationary state corresponding to the generalized Gibbs 
ensemble and next relaxes to the true steady state. This 
is what we numerically observed in Fig.l (b) and (c). 
Remarkably, in Fig. 1 (b) and (c), {(jf(t)) « —0.65 in 
the quasi-stationary state, which is close to Trpiq^ erf 
(the dotted line in Fig. 1 (b)) at the inverse tempera¬ 
ture /? = 0.85 that is determined from the expectation 
value of Hp'^ at t = 10. This fact indicates that the 
quasi-stationary state is actually described by pi^ in this 
model. 

In the case (ii), on the other hand, the relaxation pro¬ 
cess towards the quasi-stationary state and that towards 
the true steady state are indistinguishable, and hence 
stable quasi-stationary behavior is not observed in the 
initial stage of relaxation. 

Because the timescale of energy relaxation becomes 
longer exponentially as the frequency increases, we ex¬ 
pect to find the case (i) for sufficiently high-frequencies 
unless there is some special reason such as conservation 
laws [15, 17], strong quenched disorder [18, 19], diverging 
timescale due to quantum criticality [48], and so on. 

Our analysis deals with general spin models, which 
makes clear why the heating is slow in a precise manner 
and leads us to the universal scenario of relaxation pro¬ 
cesses. However, in our evaluation, the single-site energy 
is overestimated and the effect of quantum interference 
is underestimated. Hence, the divergence of the FM ex¬ 
pansion presumably begins at a higher order than our 
estimation uq « l/{8gkT). We stress that our estima¬ 
tion on the timescale is a rigorous lower bound that can 
be exponentially large in frequency, and hence the actual 
timescale of the heating will be longer than our estima¬ 
tion [49]. In order to obtain a quantitatively accurate 
estimate for a specific model, we will have to study the 
quantum dynamics of the given model numerically. 

Related to the above remark, we emphasize that our 
result does not tell us about the true steady state. It 
should be an infinite-temperature state if the Floquet 
ETH holds. However, another possibility is not excluded; 
there might be an energy-localized phase [16] with van¬ 
ishing energy-absorption rate. It is an open problem to 
understand the precise condition of the Floquet ETH. 

Summary. — In summary, we have considered the 
quantum dynamics of general driven spin systems that 
have at most /c-body interactions and a bounded single¬ 
site energy g. We have rigorously shown the Theo¬ 
rem stating that the truncated Floquet Hamiltonian is 
a quasi-conserved quantity and the rate of energy ab¬ 
sorption is exponentially small in frequency. This find- 
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ing enables us to classify the initial stage of relaxation. 
It is emphasized that we need not assume short-range in¬ 
teractions in the Hamiltonian (7). For instance, = 
6a^-yJ/N, which corresponds to the Heisenberg all-to-all 
couplings, satisfies the condition of Eq. (6) with a fixed 
value of g even in the thermodynamic limit. Therefore, 
the result in this paper is applicable to most physically- 
relevant spin models. However, as seen in Eq. (6), our 
argument excludes bosonic systems. We expect that our 
analysis will help to understand even for bosonic systems. 
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Note added. Recently, closely related results in a dif¬ 
ferent approach have appeared [36, 37]. 
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NUMERICAL DEMONSTRATION OF THE EM EXPANSION 

We consider the same model that we discussed in the main text. Namely we take a toy model where the dynamics is 
driven by H {t) = for the first half period and H (t) = for the second half period. The time-evolution operator 
for the single period U is given by 

We show the time-evolution of (cf (t)) for several system sizes and periods in Fig. SI. The initial state is chosen as 
the all spin-down state. Each graph of Fig. SI shows that (erf (t)) remains finite when the system size is small, while 
it eventually tends to zero when the system size is large. This implies that a large system eventually reaches the state 
of infinite temperature, which is consistent with the Floquet ETH. By comparing the three graphs of Fig. SI, we find 
that the timescale of the heating becomes longer as the period of the external field decreases. This is also consistent 
with our result on the timescale of the heating. 

The arrows in the figures indicate the expectation value of erf in the equilibrium state of the Oth-order truncated 
Floquet Hamiltonian , that is, (crf)qss = TrpiqVf with piq^ = /Tre^^^r” . The temperature is determined 

from the expectation value of at t = 10. In this case, fi = 0.85 and {erf )qss = —0.65. We can see that in the first 
stage of the relaxation, (erf (t)) relaxes to (erf )qss, which indicates that a quasi-stationary state observed in this model 
is actually described by the equilibrium state of the truncated Floquet Hamiltonian. 

We next demonstrate the numerical evidence that the FM-expansion can diverge in this model. In Fig. S2, we show 
the operator norms of as a function of n. It looks convergent up to some n but it starts to diverge when n 

exceeds this value. As far as we computed, the divergence of the FM expansion occurs for T > 1.0 when N < 24. 

It is emphasized that the divergence begins at a higher order term than our estimation uq = \l/(8gkT) — IJ. 
Although the theorem presented in the main text becomes meaningless when l/(8gkT) < 1 because uq becomes 
negative, the numerical simulations given here and in the main text show that exponentially slow heating is observed 
even in that case. From this observation, we must conclude that although our rigorous analysis clearly shows that the 
timescale of heating is extremely slow in the high frequency regime for general spin systems, it does not give accurate 
quantitative details. Numerical investigation is necessary in order to make a quantitatively accurate prediction for a 
given specific model. 
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FIG. SI. (color online) Time-evolution of af for (a): T = 1.02, (b): T = 1.2, and (c) T = 1.5. The arrows indicate the 
equilibrium value with respect to H^\ i.e. Trp^af with /Tre”^^^ . The temperature is determined by the 

expectation value of in the initial state {/3 = 0.85 in this case). 
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FIG. S2. (color online) The operator norm of nth order term of the FM expansion (left) and the nth order truncated Floquet 
Hamiltonian. 
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PROOF OF THE THEOREM 
Preliminary 

Here, we provide the proof of the theorem. As a preliminary, we define the “/c-locality” and the “g-extensiveness” 
as the properties of operators. We say that an operator A is fcA-local if it is decomposed as 

A= ^ ax, (S2) 

X:IXI<kA 

where A is a set of lattice sites and |A| is the number of the sites in X, see setup in the main text. This operator A 
is said to be g^-extensive if 


\Wx\\<gA 

X■.X^i,\X\<kA 


(S3) 


for all lattice sites i = 1, 2,..., A. According to this definition, the Hamiltonian considered in this paper is assumed 
to belong to the class of fc-local and g-extensive operators. 

Let A be a fc^-local and g^-extensive operator and let i? be a fc^-local and gs-extensive operator. Then we readily 
find that the commutator of A and B is {kA + A:B)-local and [2(fc^ + fcB) 5 ^ 5 B]-extensive. Using this iteratively leads 
to that the multiple commutator in Eq. (8) is (n + l)fc-local and [g{2gk)'^{n + l)!]-extensive. From this we find that 
is (n + l)fc-local and g„-extensive with 


9n 


{2gk)‘^n\ 
n + 1 


g- 


(S4) 


We can derive useful inequality for the (multiple) commutators of /cy^-local and ^yi-extensive operator A and fes-local 
operator B = E|x|<feB bx: 


\\[A,B]\\<2gAkBB, 


(S5) 


where B = J2\x\<kB ll^^ll- More generally, we can show that for any fc^^-local and -extensive operators {Ai\, 


||K,K_i,...,[Ai,H]...]]|| <i?n2g^,A,, (S6) 

2=1 


where Ki = kB+ ^Ai ■ 

By applying this inequality to Eq. (8), it is shown that the coefficient of the EM expansion can be decomposed 
as Vln = J2x-\x\<{n+i)k '^x, where wx is an operator acting on the sites in X, with 

n;r"= ^ ||n;x||T"<2gjVy ^^g^_^j"^’ . (S7) 

X-.\X\<(n+l)k ' ' 

Obviously ||H„|| < and thus Eq. (9) is also derived. 


Time evolution of local operators over single period 

In order to prove the theorem, we first consider the time evolution of a local operator O in the single period. The 
exact time evolution of the operator O in the Heisenberg picture is given by 

0{t) = U\t)OU{t) = re''^odt'L{t')Q^ 

where L{t){-) = [H{t)^-] is the Liouville operator. We define the approximate time evolution under the truncated 
Floquet Hamiltonian H^°'^ as 


Qino) ^ °^*0, 


(S9) 
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where L^p°\-) = [H^°\ •]. We can show the following lemma: 

Lemma. Assume that Hq, V{t), and H{t) are k-local and g-extensive. Then, for an arbitrary (/ + l)k-local operator 
O = J2\x\<{i+i)k^^ period T < l/{^gk), the following inequality holds: 

\\o{T) - o("“)(r)|| < (sio) 

where O = X]|x|<( 7 +i)fc = \yl{^g^T) — IJ. In particular, for O = Hq, the following stronger bound 

exists: 

\\Ho{T) - < 8g^k2-^°NvT, 

where Nv is the number of sites subjected to the periodic driving. 

Note that Eqs. (S8) and (S9) are rewritten as follows 

0{T) = Te^^o dtmo = Y,n=o T^AnO, 

O("o)(r) = ^ r”A^r^O , 

where L(t) = [H(t), ■] and = [hA°\ and An is the coefficient in the Dyson expansion given as 

dtn--- [ dtiLftn). ■ .L{ti). (S13) 

Jt2 

The coefficient A^°'^ is given by 

n .^ / r \ 

E ':^x(Y.^k + i)=n\Li,...U^, (S14) 

r-=0 ■ \i=l / 

0</i<no 

where L/ = [fl;,-] and x is the indicator function defined by x(True) = 1 and x(False) = 0. The truncated FM 
expansion is exact up to 0{T'^°) in the sense that 

An = A^°'^ for all n < uq. (S15) 



(Sll) 


(S12) 


By using this fact, we have 

OO 

\\o{T)-d{T)\\< ^ r"(M„o|| + p^io||). 

n—nQ-\-l 


(S16) 


We bound ||.T„0|| and ||.A^r“^0|| from above, using Fq. (S6). First, we show from Fqs. (S13) and (S6) that ||.4„0|| 
is bounded as 

||A.O||<- sup \\L{tn)...L{h)0\\<{2gkr^-^^0. (S17) 

By using (n + /)!/[n!/!] < 2"+^ and AgkT < 1/2, we have 

OO 

r"||^„0|| < 8g/cT2^-"'>a (S18) 

n—no-\-l 

For O = Ho, we can give a stronger bound. Because L{ti)Ho = [H{ti),Ho\ = [F(ti),i7o] = —iol^(ti), where 
Lq = [Ho, •], we have 

\\AnHo\\ < ^ sup ||L(t„)...L(t 2 )LoF(ti)||. 
n'. o<ti...,t„<T 


(S19) 
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Because V{t) is fc-local, it is written as V{t) = X]|js:|<fc where vx{t) is an operator acting non-trivially only on 

the domain X. By almost the same calculation, we have 

OO 

^ T^\\AnH4<AgkT2—°Vo, (S20) 

n=no + l 


where Vq = snpg<«y ||nx(t)||- Because V{t) is g-extensive, Vq < gNy if the driving is applied to Ny sites. 

Next, we evaluate an upper bound of 


E -,x{Y.(h + l) = n]\\Li,...LiM- 


^=0 {iJLi 

0</i <no 


(S21) 




(S22) 


By using Eqs. (S6) and (S4) with gi < {2gknoyg, because l\ < Uq for I < no, we obtain 

\\Li, ... LiM < {2gkn2gknor-^j^-^-^0. 

[n + I — r)l 

By using (n + /)!/[r!(n + I — r)!] < 2"+^ and ^ + i) = n) = {n — l)!/[(n — r)!(r — 1)!], we obtain 

\\A(^o)o\\ < 2-+^ (2gfc)’'(25A:no)"-’'0 = 2^[4gk{no + l^-^^gkO. (S23) 

Because no is the maximum integer not exceeding 1/iSgkT) — 1, we have AgkT{nQ + 1) < 1/2 and thus 


< 8gkT2^-'^°0. 

n—no-\-l 


(S24) 


By substituting Eqs. (S18) and (S24) into Eq. (S16), we complete the proof of Lemma (SIO). 

Eor O = Hq, we have Li^Hq = —LqXIi^ and thus \\Li^ ... Li^HqW = ||L/^ ... Li^_^LoXli^\\. Clearly, Ir = 0 has no 
contribution, and for 4 > 1, we can use Eq. (S7). As 1 < < no, we have < {2gknoy^gNy/2. Applying Eq. (S6) 
and the upper bound of we obtain 


Y T"||A/"“)iLo|| < 4g2/cTAfy2-"«. 

n—nQ-\-l 

Using Eqs. (S20) and (S25), we obtain Lemma (Sll). 


(S25) 


Proof of Theorem 

Erom the relation HY\t) = hY\ "we note the inequality for any positive integer m: 

\\HY\mT) - hY\\ < m\\HY\T) - H^'^W 

no no 

< m^T"||L!„(T) - Ui"“)(T)|| = m||no(r) - no(T)|| + m^ T"||L!„(r) - Ui"«)(T)||. (S26) 

n—0 n—1 

Applying Lemma and using Eq. (S7), we obtain 

||no(r) - HoiT)\\ < 8g^k2-^°NyT, (S27) 

||0„(T) - nyy°'>{T)\\T^ < 16gkTTy2-<'^°-^'>T^+^ < 32g^k2-^°NyT . (S28) 

(n + 1)^ 
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By using n! < tiq < l/{8gkT)^ and l/(n + 1)^ < 1/4 for 1 < n < uq, we have 

no oo / ^ \ n 

^ ||ll„(r) - Oi"'>)(T)||T" < 8g^k2-^«NvTj2 U ) 

= 8g^k2-'^«NvT. (S29) 

By combining this and Eq. (S27) with Eq. (S26), we obtain the theorem 

\\H^p°\t) - < 16g^k2-^°Nvt, (S30) 

where t = mT. This completes the proof of Theorem. 


